Enhancing the intense field control of molecular fragmentation 
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We describe a pump-probe scheme with which the spatial asymmetry of dissociating molecular 
fragments — as controlled by the carrier-envelope phase of an intense few-cycle laser pulse — can 
be enhanced by an order of magnitude or more. We illustrate the scheme using extensive, full- 
dimensional calculations for dissociation of and include the averaging necessary for comparison 
with experiment. 



In recent years, considerable experimental effort has 
been invested in developing the ability to control chem- 
ical reactions with intense, few-cycle laser pulses [1-3]. 
The canonical reaction chosen has been molecular dis- 
sociation, and a common measure of the degree of its 
coherent control is the spatial asymmetry of fragments 
relative to a linearly polarized laser field. Quantum me- 
chanically, this asymmetry arises from the interference of 
pathways that lead to even and odd parity states [4-6]. In 
strong fields, these pathways can involve many photons, 
and the relative phase between pathways — and thus the 
outcome — can be controlled by varying laser parame- 
ters such as the carrier envelope phase (CEP) [1-3, 5-14] 
or the relative phase between different colors [15-22]. In 
this Letter, we will focus on control via the CEP. 

In the dipole approximation, the CEP Lp for a Gaussian 
laser pulse £{t) is defined from [23] 

£(t) = Soe-^^^^' cos{ujt + if). (1) 

Generally, the largest CEP-dependent asymmetries have 
been observed for ionized electrons [1]. The asymme- 
tries for the nuclear fragments resulting from dissociation 
have, unfortunately, been much smaller [3, 8, 9]. These 
weak effects — combined with the ongoing challenge of 
producing intense, few-cycle, CEP stabilized pulses — 
greatly limit experimentalists' abilities to measure and 
explore this intriguing means of control. One important 
recent advance is the ability to measure the CEP of each 
pulse [24, 25], alleviating the need for CEP stability dur- 
ing the measurements. 

While CEP-dependent asymmetric break up of II2" was 
predicted a few years ago [2], successful measurements 
have not yet been made starting directly from this bench- 
mark system, e.g. in an ion beam experiment [26]. Ex- 
periments have instead begun with the more complicated 
H2 [3, 9]. With only one electron, the number of control 
pathways for H2" is smaller than for H2 making the inter- 
pretation more straightforward. Moreover, the theory at 
sub-ionization intensities can be done essentially exactly. 

The technical challenges of an ion beam experi- 
ment [27] are obvious reasons that the II2" experiments 
have not yet been done. A more fundamental problem, 
however, lies in the fact that H2" typically comes in a 
broad rovibrational distribution in such experiments [27]. 
Unfortunately, dissociation of II2" from different initial v 



by a linearly polarized laser pulse gives fragments with 
similar energies. Since the asymmetry produced by each 

V is slightly different, the incoherent averaging over ini- 
tial V required for an II2" beam tends to wash out the 
overall asymmetry [28]. Moreover, one-photon dissocia- 
tion of higher v dominates the total signal. And, since 
one-photon dissociation produces a single nuclear parity, 
its momentum distribution is symmetric, masking the 
desired asymmetry. 

In this Letter, we present a scheme to greatly enhance 
CEP effects. This enhancement is largely achieved by de- 
pleting the undesired higher-^ states with a long, weak 
pump pulse. Subsequent dissociation of this prepared 
state by a few-cycle probe pulse gives a momentum distri- 
bution with an order of magnitude enhanced asymmetry 
compared to that of an initial incoherent Franck-Condon 
distribution. In fact, our scheme gives larger asymmetries 
— at longer pulse lengths — than have been observed so 
far in II2 experiments [3, 8, 9, 13]. We also propose ways 
to separate the pump and probe signals. To support our 
claims, we present theoretical CEP-dependent p+H mo- 
mentum distributions in addition to the up-down asym- 
metry. Calculating such a differential observable — along 
with averaging over the intensity distribution of the laser 
focus — permits us to quantitatively predict the experi- 
mental outcome and to provide deeper physical insight. 

To obtain the momentum distribution, it is neces- 
sary to account for the nuclear rotation. We thus solve 
the time-dependent Schrodinger equation in the Born- 
Oppenheimer representation, including nuclear rotation, 
nuclear vibration, and electronic excitation, but neglect- 
ing ionization as well as the Coriolis and all nonadiabatic 
couplings. The nuclear rotation is included as an expan- 
sion of the wave function over the total orbital angular 
momentum (J) basis (see [29] for details). 

To prepare the system, we use a 785 nm, 45 fs long 
pump pulse with an intensity of 10^^ W/cm^. This 
relatively long, weak pump pulse depletes the higher v 
states, eliminating their spatially symmetric dissociation 
signal. Figure 1 shows the ^-distribution before and af- 
ter this pump pulse. In an incoherent Franck-Condon v- 
distribution, 9.58% of the population lies in 'u > 8. Since 
90% of this population dissociates in the pump pulse, 

V > 8 becomes only 1.36% of the total remaining bound 
population, ensuring that their contribution to the dis- 
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FIG. 1: Vibrational state distribution before and after the 
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45 fs, 10 W/cm pump pulse. 



sociation signal by any subsequent probe pulse will be 
negligible. Consequently, we performed probe pulse cal- 
culations only including v =0-7, after verifying for a rep- 
resentative case that v > 8 affected the asymmetry by 
much less than 1%. 

To quantify the enhancement, we compare the re- 
sults from an initial incoherent Franck- Condon distri- 
bution interacting with only the probe pulse ("probe- 
only" ) to the signal from the probe part of our proposed 
pump-probe scheme ("pump-probe"). We used a 7 fs, 
785 nm probe pulse in both cases. In the pump-probe 
scheme, all calculations were performed at a fixed time 
delay of 267 fs unless stated otherwise. Since ioniza- 
tion is neglected, we limit the peak intensity to no more 
than 1.2 X 10^^ W/cm^. For the peak intensities above 
10^^ W/cm^ required for intensity averaging, our calcu- 
lations included p+H(2/) manifold in addition to Isag 
and 2pau channels. The total population of the p+H(2/) 
states was less than 5% even for the highest intensity. 
Consequently, we present momentum distributions based 
on just the Isag and 2p(Ju channels. 

The fundamental physical observable we focus on is 
the p+H relative momentum distribution p(K), which is 
the most differential observable in recent experiments in- 
volving H2" dissociation [26, 27, 30]. To calculate p(K), 
we project the final wave function onto scattering states 
that behave as exp(iK • R)^is^ asymptotically, where 
R points from proton A to proton B and (j)isA is the 
hydrogen ground state wave function centered on pro- 
ton A. The momentum K thus points from H to p. 
This scattering state, with the nuclear spin included, is 
then symmetrized to account for the identical nuclei [31- 
33]. Finally, the momentum distribution [or its energy- 
normalized equivalent p{E^K) with E = i^^/2/i, /i the 
nuclear reduced mass, and K = {Ok-,^k) the direction 
of K with respect to the polarization direction] is 
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FIG. 2: Franck- Condon- averaged K^p(K) (integrated over 
(Pk) for probe-only [p(K) is reflected to —K± for clarity] for 
(a) = 0, (b) = 7r/4, (c) ip = 7r/2, and (d) p = 37r/4 ( Gray 
dotted lines mark K = 7.5, 10, 12.5, and 15 a.u.). (e)-(h) 
are same as (a)-(d) but for pump-probe. All cases used a 7 fs, 
10^^ W/cm^ probe pulse. 



with 



Cjp = {-iye-''^'-{EJp\Fjp{tj)), p = g, 



u. 



(3) 



Here, \Fjp{tf)) are the Isag and 2pau nuclear radial wave 
functions at the final time tj, while \EJp) and 5jp are the 
corresponding energy-normalized scattering states and 
phase shifts, respectively. 

Equation (2) shows that although a linear combina- 
tion of Isag and 2pau is necessary to localize the electron 
as an atomic rather than a molecular state, the spatial 
asymmetry of p+H is due to the interference of even and 
odd nuclear parity states. This distinction is brought 
into sharp relief when nuclear rotation is included in the 
calculation since using simply lsag±.2pau nuclear wave 
function — as is done in calculations without rotation 
— would produce two distinct p+H momentum distri- 
butions, where clearly only one can be measured. It is 
the symmetrization requirement that dictates the proper 
coherent combination to use. This issue is not new, how- 
ever, and always arises for identical particle scattering 
where it is known that the primary differences occur for 
^x~7r/2. Since intense-field dissociation of H2" produces 
very few fragments at this Ok^ the consequences of ana- 
lyzing incorrectly are less pronounced. For more compli- 
cated systems, however, this need no longer be true. 

The Franck- Condon- aver aged momentum distribu- 
tions for several CEPs are shown in Figs. 2(a)-(d) for 
the probe-only case and in Figs. 2(e)-(h) for the pump- 
probe case. The momentum distributions in all cases ex- 
hibit preferential alignment along the laser polarization. 
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FIG. 3: (a) Asymmetry defined in Eq. (5) for tiie probe-only 
case and (b) for tiie pump-probe case for a tfwhm = 7 fs and 
/ = 10^^ W/cm^ pulse. 

Moreover, since tfie energy distributions 

p{E) = j p{E,K)dnK (4) 

for the Isag and 2p(Ju channels of individual vibrational 
states overlap roughly in the range 0.5-1.5 eV, we expect 
spatial asymmetries to appear roughly for 6<K<10 a.u. 
The results shown in Fig. 2 for both experimental scenar- 
ios are consistent with this expectation. The momentum 
distribution for = + tt is the mirror image of the 
momentum distribution for cp^ as guaranteed by the fact 
that cos{ujt + tt) = — cos ujt in Eq. (1). 

While the two experimental scenarios clearly show 
qualitative differences, the strikingly different distribu- 
tions make it difficult to judge which produces the larger 
asymmetry. We thus turn to the quantitative measure 
of the asymmetry used in previous studies [3, 6, 8]: the 
normalized asymmetry parameter A{E^ip), 

A{E, if) = piE)-' [p(£;)up - /9(i?)Dow„] . (5) 

For simplicity, we integrate over the whole upper and 
lower hemispheres in the "Up" and "Down" distribu- 
tions, respectively, although a narrow angular cut along 
the laser polarization direction might be chosen to en- 
hance A as in some experimental studies [3, 8]. Figure 2 
shows why such cuts are effective since the strongest CEP 
dependence lies at small Ok- Although the total energy 
spectrum p{E) in principle also depends on CEP [6], we 
found negligible CEP-dependence in the Franck-Condon 
averaged p{E) and thus expect essentially no contribu- 
tion to the CEP dependence from the denominator of 
A. 

Figures 3(a) and 3(b) show A{E^ (p) for the proble-only 
and the pump-probe, respectively. For the probe-only in 
Fig. 3(a), we can already see reasonable asymmetry in 
the range 0.2-2.5 eV where it oscillates between —0.12 
and 0.12. Comparing Fig. 3(a) and Fig. 3(b), however, 
we find a five- fold enhancement of \A{E^ (p) \ in the pump- 
probe case for this intensity. 

The most crucial factor determining whether an 
intensity-dependent effect is experimentally observable 
is whether it survives intensity (or focal volume) av- 
eraging. We thus intensity-averaged our results for a 



7 fs laser pulse with cp = and a peak intensity of 
1.2 X 10^^ W/cm^. We used the two-dimensional ge- 
ometry of Ref. [30] to perform the intensity-averaging 
following the procedure described in [28]. Figures 2(a), 
2(e) and 3, show a clear up-down asymmetry for = 
in both probe-only and pump-probe cases, and we will 
check if it survives intensity averaging. Figure 4(a) shows 
A{E, 0) before and after the intensity averaging for both 
cases. The corresponding total KER distributions are 
plotted in Fig. 4(b) to show that the small p{E) is the 
reason for the large A{E^ 0) at higher energies. In the 
pump-probe case, the intensity averaging has only been 
performed over the probe-pulse intensity distribution as- 
suming that the weak pump intensity can be made uni- 
form across the probe focal volume. 

For the probe-only case, intensity averaging reduces 
A by more than a factor of three over the entire energy 
range shown in Fig. 4(a) and makes it 17 times smaller for 
0.5 to 1.0 eV, where p{E) is large and the single-intensity 
A is largest. This significant reduction in A is due to 
the fact that one-photon dissociation — which shows no 
asymmetry — can occur for v > 7 at very low intensi- 
ties (^10^^ W/cm^). These symmetric contributions are 
thus amplified by the intensity averaging and swamp any 
asymmetry because essentially all v contribute to these 
KER. Figure 4(a) thus shows that intensity averaging 
makes it very challenging to measure CEP-effects for a 
single 7 fs or longer pulse in an experiment. 

For the pump-probe case, A is also reduced from the 
single intensity value — but to a much lesser extent than 
in the probe-only case. In fact. Fig. 4(a) shows that 
even after intensity averaging A is an order of magni- 
tude larger using the pump-probe scheme compared to 
the probe-only results. Moreover, we have found that the 
pump-probe scheme produces a CEP-dependent asym- 
metry after intensity averaging even for 10 fs pulses. 
These pulses are much longer than the 6 fs pulses that 
have been used to date to observe CEP effects [3, 9]. 
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FIG. 4: (a) Asymmetry from the intensity- averaged p(E, K) 
for the probe-only (thin dashed lines) and the pump-probe 
(thick dashed lines) cases as well as for a single 7 fs probe 
pulse with a peak intensity of 1.2 x 10^^ W/cm^ (thin and 
thick solid lines, respectively), (b) KER distributions for the 
cases shown in (a) normalized to the same peak value. 
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Besides depleting the high- lying vibrational states, 
the pump also impulsively aligns the molecule [34- 
36]). To investigate the sensitivity of A to the align- 
ment, we calculated the asymmetry for three different 
pump-probe time delays with aligned (cos^ ^)=0.56, anti- 
aligned (cos^ 6>)=0.22, and dephased (cos^ 0)=0A0 angu- 
lar distributions ((cos^ 6>)=l/3 for an isotropic distribu- 
tion) . We found that the maximum A was largest for the 
aligned distribution, followed by the ant i- aligned, with 
the dephased smallest. The enhancement of the aligned 
A over the dephased was roughly 30%. For this reason 
we have shown here calculations for the 267 fs delay cor- 
responding to the aligned distribution. This exercise also 
served to establish that the major source of the ten-fold 
CEP-dependent asymmetry enhancement is the deple- 
tion of the higher-'u states. 

Another concern for experimentally observing the pre- 
dicted enhancement is the fact that in our pump-probe 
scheme the pump pulse already produces fragments. So, 
separating the probe signal from the pump signal is cru- 
cial. Although the dissociating fragments from both 
pulses overlap in momentum, we expect the asymme- 
try will still be large in a combined pump-probe sig- 
nal for two reasons. First, the momentum distribution 
from the long pump pulse exhibits narrow peaks corre- 
sponding to higher vibrational states. Therefore, in the 
combined pump-probe momentum distribution, the sym- 
metric structure would be very localized in KER, giving 
small overlap with the broad asymmetric signal resulting 
in larger asymmetry than the probe-only case. Second, 
we found that preparing the initial state greatly increased 
the dissociation probability of the lower vibrational states 
for aligned (2.33 fold) and dephased (1.74 fold) pump- 
probe cases compared to the probe-only, thereby enhanc- 
ing the ratio of the asymmetric signal to the symmetric 
signal. 

The contrast between pump and probe signals can be 
further improved over the present case using pump pulses 
longer than 45 fs, thus increasing the depletion of the 
higher vibrational states and making the pump signal 
even more structured. A longer pulse will give more 
alignment, which might also enhance asymmetry. Ad- 
ditionally, instead of using the whole upper and lower 
hemispheres to define A^ an angular cut can be used to 
isolate the aligned asymmetric distribution. 

As the dissociating fragments primarily lie along the 
laser polarization, it might be better to use orthogo- 
nal laser polarization directions for the pump and probe 
pulses to separate their signals [13]. For this, one might 
want to use the time-delay when the molecules are an- 
tialigned relative to the pump polarization to improve the 
signal. A circularly polarized pump pulse could also be 
used. Since depletion is the major reason for enhanced 
CEP effects, we believe the effect will survive using differ- 
ent laser polarizations. An intensity differencing scheme 
might also be useful to enhance asymmetry [37]. 



In this Letter, we have presented a prescription for sub- 
stantially enhancing CEP effects on the spatial asymme- 
try of intense-field-induced fragmentation. We have il- 
lustrated our proposal with essentially exact calculations 
for the benchmark system H2", including the important 
averaging implicit in experiments and found a ten-fold 
increase in the asymmetry. In addition, we have sug- 
gested several other steps that could further increase the 
asymmetry. 

While even greater enhancement could be realized by 
preparing the system in a single initial our scheme pro- 
vides a more easily followed experimental avenue yield- 
ing a narrow vibrational state distribution. We believe 
that our scheme is equally applicable to neutral molecules 
with studies focussed on dissociative ionization channels, 
where the initial intense long pump pulses can serve to 
ionize, dissociate undesired vibrational states, and align 
the molecular ions. 
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